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Abstract 


In this thesis we have implemented some well known techniques for solving queueing 
networks that have only non blocking nodes in the network m a user friendly software 
package called Queueing Network Analysis Tool (QNAT) We have also developed some 
techniques to reduce fork join queues that may or may not have synchronising queues and 
have made them available in QNAT QNAT can solve open and closed queueing 
networks with single class jobs and open closed and mixed networks with multiple class 
jobs QNAT does not allow jobs to change classes 

For the analysis of open networks we use the technique developed for networks of 
GI/G/m queues by Whitt For closed networks we use the well known Mean Value 
Analysis (MV A) technique For mixed networks networks that have open and closed 
classes of customers we first find the utilisation of the servers for the customers of the 
open classes This utilisation is used to reduce the service rate (inflate the service 
demand) for the closed classes of customers We now consider only the closed classes of 
customers but with the reduced rate of service at the queues and then solve the network 
using the MVA technique This gives us the performance parameters for the closed 
classes of customers The performance of the open classes of customers is then obtained 
from those of the closed classes of customers 

In QNAT fork join queues can also be specified in the network The fork join queue may 
or may not have synchronising queues after the service queue in the siblings The network 
may be open or closed In all these cases a flow equivalent service centre (FESC) for the 
fork join queue is obtained using suitable techniques The fork join queue is replaced by 
this FESC and the resultant network is solved using standard techniques Fork join 
queues with synchronising queues is not supported by QNAT for open networks QNAT 
can also handle queueing networks that have finite capacity queues and its 
implementation are detailed in a companion thesis 



QNAT has a user friendly graphical user interface (GUT) based front end from which we 
can create a network save the network with suitable names 0|>en the networks that were 
previously saved print the solution results into a file etc QNAT also has an online help 
to guide the user at every stage The executable files have been developed in such a way 
that future modifications like providing a simulation option or adding any new 
capabilities can be easily incorporated 

A beta version of QNAT is being distributed over the Internet It can also be obtained 
from the ERNet/Telematics Laboratory of ITT Kanpur 
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Chapter 1 


Introduction 


The design and development of many of the rapidly evolving complex systems requires 
that their throughput and delay performance be studied during the design and 
development process Examples of such systems include computer systems 
manufacturing systems networked communication systems and distributed datibise 
systems There are many ways to study the performance of such systems An obvious 
method would be to use intuition and perform trend extrapolation The use of such i 
method relies heavily on the experience of the performance analyst and the reliibility 
associated with this is very low because of the absence of a scientific basis for the 
conclusions from the study Further there will always be very few people who c in 
develop the ability and the intuition to accurately predict the performance ot systems 
through trend extrapolation A second option available to performance analysts is th it ol 
experimental evaluation Here the actual system will be subjected to vinous loul 
conditions and the relevant performance parameters measured This appro ich is 
obviously very expensive Further it might not be possible to obtain realistic conditions 
for the experimental setup A considerable amount of extrapolation may be neccss irv 
from the results of a few tests and the results obtained by this approach may not be easily 
generalised A third approach the one preferred by many performance analysts is th it ot 
analytical (or mathematical) modelling Among the many mathematical techniques th it 
can be used to model systems let us consider queueing network based models 

Many systems can be represented as a network of queues Jobs enter a queue wait foi 
their turn to receive service at the queue and upon completion of service they may either 
enter another queue or leave the system The arrival process to a queue the service 
discipline at a queue the service process at a queue and the routing probabilities ittcr 



service at a queue are some of the parameters that specify the queueing network mode! of 
a system Under certain assumptions of the parameters specifying the queueing network 
an analytical solution may be obtained by evaluating a set of linear equations obtained 
from the queueing network specification If the analytical solution is not easily available 
a computer simulation model for the queueing network model of the system untier stud) 
can be created The principal advantage of a simulation study is its flexibility However 
simulation models are generally expensive to obtain because of the need to write and 
debug complex computer programs 

From the above discussion it is evident that analytical modelling is the preferred method 
in the performance analysis of systems In this thesis our interest is in queueing network 
based analytical modelling of systems 

Queueing network models with simple assumptions achieve a fivourable bal inee 
between accuracy and efficiency These models achieve relatively high accuriey at i 
relatively low cost queueing network modelling is a small subset ol the tcthniqiius ol 
queueing theory Much of queueing theory is oriented towards modelling a coinplex 
system using a single queue or a service centre with complex characteristics A detailed 
performance analysis of such systems can be obtained by using sophistic Ued anilyticil 
techniques Simulations can also be performed on such complex single qiiciii 
representations of systems but the performance results obtained are based on complex 
assumptions and will generally be difficult to generalise In queueing network modelling 
a network of simple queues are considered and simple evaluation algorithms are used to 
obtain meaningful performance measures with an appropriate balance between accuricy 
and efficiency 

Constructing a queueing network model involves the following steps delinition ol tlu 
system parameterization and evaluation Defining systems is relatively easy whcncvti 
these systems can be represented as a network of queues In most systems that interest us 
parameterization of the input variables that define the system and defining the output 



variables of interest is relatively straightforward and system description is e isily 
accomplished It is the solution of the model that requires extensive m ithem ilic il 
calculations Further if the system is to be evaluated for a number of different sets of 
input parameters the system evaluation becomes tedious and time consuming Thus 
there is a great need for a tool to solve queueing network and obtain their perform incc 
parameters 

The steps in the specification and solution of the queueing network model for any system 
can be integrated into a single user friendly package where the simple terminology of the 
queueing network system rather than the language of queueing theors is used Wc au 
aware of the existence of a large number of software packages that base been dc\c loped 
to specify and solve queueing network systems Some of these are CADS (1978) BfST/l 
(1980) CAPTURE/MVS (1982) RESQ (1982) MAP (1982) PANACEA (1982) 
CMF/Model (1983) and QNA (1983) ^ ^ Most of these packages are not easil\ isail ihle 
Also in these packages the system specification is through a text file which in these 
davs of GUI based user friendly computing would not be very useful Further all ihesi 
packages have been designed to solve only a restricted class of queueing networks In 
1996 an integrated user friendly package called NAST Network Anahsis and Simulation 
Tool was developed in the ERNet/Telematics Laboratory of ITT Kanpur'’’ This is 
similar in its capabilities to the packages mentioned above except that the input s\siem 
and Its parameters is specified through a GUI and the softwaie is de\ eloped loi Mierosoli 
Windows™ on PCs and uses Mathematica™ as the compute engine In this thesis we 
ha\e developed an updated version of NAST which can sohe mans new queueine 
models like BCMP queueing networks and queueing networks with tork join queues in 
addition to those supported by NAST This package is called QNAT (Queueing Network 
Analysis Tool) QNAT can also solve queueing networks that ha\e finite capacity nodes 


’ Edward D Lazowska JohnZahorjan G 
Performance Prenliee Hall Inc 1984 


Seott Graham and Kenneth C Seeeik Quantitciti\i Swum 


^ ard Whitt The Queueing Network Aiicilwer BSTJ Vol 62 No 9 Nov 198 
’M N Umesh A Window based Package for Queueing Network Aimlwii MTeehThcMs IlTKuipur 
Feb 1997 



The following briefly summarises the contents of this thesis In Chapter 2 we give a brief 
overview of the queueing network models that we have incorporated into QNAT In 
Chapter 3 we provide the details of the algorithms used in this package We will not 
discuss the algorithms for solving queueing networks that have finite c ip icily nodes 
These will be discussed in a companion thesis’ We also provide examples and 
verifications for the algorithms Chapter 4 contains an overview of the contents of the 
QNAT software package conclusion and outline for future work 


‘ H Tahilramani Sohmg Queueing Networks nith FintU Capcicin NacUi m QNAT MTtch Thesis In 
preparation 
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Chapter 2 


Queueing Network Models 
2 1 Introduction 

Most complex systems that provide any kind of service at multiple nodes can be modeled 
as a network of queues by an appropriate definition for the arrival process the servat 
process and the routing of the customers after service at a node Thus by evalu iting ihs 
performance of the equivalent network of queues representation of the systsm wc cm 
study the performance of the original complex system Consider a machine shop ex impk 
with a number of machines Jobs will enter the machine shop from outside The\ will 
receive service from the various machines in the shop in some possibK random 
sequence and finall> leave the shop In the equivalent queueing network model the 
machines are represented by the individual queues providing service the jobs rcprescnisd 
bv the customers and the sequence of machines used by jobs corresponds to the routin^ 
process Under fairly general assumptions for the arrival service and routin_ processus 
the equivalent network of queues can be shown to be a multi dimensional birth dc ith 
process 

The first useful result available for queueing networks is Jackson s Theorem This 
theorem states that if the external arrivals into the queues form a Poisson process the 
service time at each of the queues is independent of the service at the other nodes and has 
an exponential distribution and if the routing is Markovian i e the probabilitv of goin^ 
from node i to node j after completion of service at node / is independent of the past then 
the joint distribution of the state vector of the system n=[ni ot n -5 Hn] has a product fomi 
solution of the form 

Pi(ni) P (n-.) PnCon) 

P(ni ns nN)= for/2=0 1 2 (2 1) 

G 


s 



where P(ni nj Hn) is the probability that the network of N queues is in state n = 
(ni n 2 nw) n, is the number of customers in queue t P,(n,) is the probability of 
queue i having n, customers and is calculated by treating node / independent of the other 
queues in the network and by obtaining the arrival rate into that queue from the arrival 
process into the network and the routing process G is the normalising const int to m ike 
the right hand side of the above equation a probability mass function 

The following example illustrates Jackson s theorem Consider the simple network ol two 
queues shown in figure 2 1 Jobs enter queue 1 with rate Xi After receiving service it 
queue 1 the customer comes back to node 1 for another service with probability p or gous 
to queue 2 with probability (1 p) Customers can also arrive to queue 2 with rite X. The 
service rates at nodes 1 and 2 are pi and p respectively 



X| 

— ► 


|X 

' — ► 

i 

L 

/, arrival rate it node i 


— ► (I) 

-T— 

TTp) (2) 


p, serv ice rale U node i 



ii 


P 



p 





Figure2 1 An 

open queueing network 



From Jackson s theorem 

P(ni n ) = Pi(ni) P (n ) = pi"' p " / G if pi < 1 and p < 1 

where 

Xi X| + X 

Pi= P-= 

(Ip)pi P 

It can also be easilv seen that the normalising constant G will be obtained from tlii. 
normalising constants of the two M/M/1 queues (1 pO le G' = (1 pOfl p ) 







Statistics of the waiting times and the queue lengths at the queues can be easily obtained 
from the above equations 

It has also been shown that if the network of queues contains a fixed number of customers 
at all times i e the network is a closed queueing network with no external arrivals into 
the network or departures from it the product form solution of Eqn 2 1 is retained In 
calculating P(n,) we have to consider a relative arrival rate rather than the absolute irriv il 
rates into the node This is because the network is closed an absolute arm il rite cannot 
be defined Also the calculation of the normalisation constant G becomes a lot mort 
difficult and was the main computational bottleneck in the useful application of this 
queueing model 

Consider the example shown in figure 2 2 There are M customers in the nctwoik A 
customer after receiving service at queue 1 will either return to queue 1 tor another 
serMce with probability p or go to queue 2 with probability 1 p After recening service at 
queue 2 the customer goes to queue 1 with probability 1 


( 1 ) 


TT^i (2) 




M Lustomers 


p serv ice rate it queue / 


Figure2 2 A closed queueing network 


The probabilit> that the system is in state n= [ni ni] is given bv 


P(ni n ) = 


Pi(ni) P 2 (ni) 
G 


= 0 


« "I ” 

pi P 


for n\+ lit = M 
for/?i+ii 


(2 2 ) 
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where Pi = Xi/^i = 1/((I p)p,) 
p2 = X / P2 = 1 / M-i 

Here we have calculated the arrival rates relative to node 2 The normalising constant G 
IS obtained by summing P(n) over all possible values 

G = 5 ; P(n) (2 

nl+nl M 

As can be seen from the above equation evaluating G requires the summation of iIk 
product terms over 'Cn terms Further these terms have to be searched m a region 
of terms This is clearly computationally inefficient and efficient algorithms arc 
necessary for closed queueing network models to be usefully applied Also as we will 
now show queue length and waiting statistics of a network of N queues with M 
customers can be expressed as functions of the normalising constant of netwoiks with V 
or fewer queues and M or fewer customers 

Let X, be the relatise arrival rate to node i |i, the service rate at node / and the X, = Xjyi 
X, is obtained by solving the equation XP=X where P the routing matrix containing thi 
probabilities P,j the probability of going to node j after service at node i G(\' M) is itu 
normalising constant evaluated for a closed queueing network with N nodes ind M 
customers The following can be easily shown 

Ltilisation at node / = p, = X, G(A^ 1)/G( A) P(n, > 1) (2 4) 

Nett arrival rate to node i = X, = p, p, (2 ‘s) 

N 

Average number of customers at node i = E(n,) = ^ (X,)" G(N 1 ) / G( V) (2 6) 

=1 

Thus we see that the solution of a queueing network model boils down to the ev ilii ition 
of a normalizing constant In the following we describe the convolution algorithm to 
evaluate G(NM) 


P(n, n, nN) = (l/G(A0)n (2 7) 

The normalizing constant can be evaluated by following the recursive algorithm 


K 



Initialize for m = 1 2 M G(0 m) =1 0 
for n = 1 2 

for m = 1 2 M 

G(n m) = G(n m 1) + Xn,G(n 1 m) 


If the relative utilization s are larger than the actual utilization s the normalizing constant 
IS greater than 1 and may be very large On the other hand if the actual utilization s are 
much larger than the relative utilization s the normalizing constant may become verv 
small This can cause numerical overflow or underflow during the calculation of the 
normalizing constant To overcome this numerical difficulty an iterative technique 
known as Mean Value Analysis (MVA) has been developed This technique avoids the 
problem of numerical overflow or underflow by dealing with the ratios of the normalizing 
constant Also it has the advantage of directlv computing the average performance 
measures The algorithm in this analysis is based on two fundamental theorems viz 
Little s Theorem and Mean Value Theorem From Little s theorem 
N.= TW.X, 

where Ni is the mean number of jobs at queue / T is the mean throughput of the network 
(on a reference queue) W, is the mean delay at queue i and X,, is the relative arrival rate at 
queue i The Mean Value Theorem^'^ states that a customer arriving to a queue in a 
product form queueing network sees precisely the same distribution of customers as an 
outside observer would see if there were one less customer in the network From these 
two theorems a recursive algorithm is developed in which waiting time of a new 
customer is computed from the waiting times computed using one less customer in the 
network Thus starting with zero customers in the network we can iterate upto M The 
outline of the algonthm is given below and is discussed in detail in the next chapter In 
this algonthm x, is the mean service time at queue i 


'M Reiser and S S La\enberg Mean Value Analysis of Closed Multichain Queueing Networks J of 
ACM Vol 27 No 2 pp 313 322 April 1980 
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Initialize for j = 1 2 1,(0) = 1 0 

for m = 1 2 M 
for i= 1 2 

= T, * /, (m 1 ) 

7’ = "I / X ^ 

=1 

for / = 1 2 ^ 

l,(m) = T*w, 


2 2 Queueing Networks with Multiple Classes of Customers 

In the previous section we discussed networks with a single class of customers Fuithcr 
we limited ourselves to FCFS queues with exponential service time distributions It was 
first shown that a queueing network model for a central server system with a processor 
sharing CPU and exponentially distributed I/O service time had product form solution il 
the CPU service time distribution had a rational Laplace transform It was also shown th ii 
the steady state behaviour of such models depends only on the mean CPU service time 
and not on its higher moments These results were then extended to arbitrary networks bv 
defining and using the concept of local balance for a queueing network to show thii 
certain other queueing disciplines like processor sharing and last come iirst served with 
preemptive resume also satisfy local balance and hence have a product lorm soluiioii 
Baskett Chandv Muntz and Palacios extended these results to a large cl iss ol 
networks that serve multiple classes of customers They showed that the product form 
solution exists for all kinds of networks open closed or mixed if thev satisfy the 
following conditions 

• All the service centres have an> one of the following four service discipline 
first come first served (FCFS) last come first served with preemptive resume 
(LCFS PR) processor sharing (PS) and an infinite server queue (IS) 

• While waiting or receiving service at a queue the customers are not allowed 
to change classes After service completion the> may change classes and 


'F Baskett K M Chandy R R Muntz and F G Palatios Open Closed and Mtxid ninunh of 
queues u ith different class of customers J of ACM Voi 22 No 2 pp 248 260 April 1975 
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should be routed according to fixed probabilities Routing is Markovian and 
does not depend on the previous visits of the customer 

• At a FCFS queue the service time distribution must be identical and 
exponential for all classes of jobs At other types of queues the service times 
should have probability distnbutions with rational Laplace transforms and 
different classes of jobs may have different distributions 

• At a FCFS queue the service time can depend only on the total queue length 
at that queue At LCFS PR PS and IS queues the service time for a given 
class of customers can depend only on the total queue length at that queue of 
that class but not on the queue length of other classes 

• In an open network the interamval time should be exponentially distributed 
Bulk arrivals are not permitted The arrival rate can be state dependtiu 

• A network can be open with respect to some classes and closed with respect to 
other classes of customers Such networks are called mixed networks 

Queueing networks satisfying the above criteria are called as BCMP networks It has also 
been shown that a queuing network that satisfies the M => M property has a pioducl 
form solution A queue is said to have M => M property if a Poisson input process results 
in a Poisson output process 

We now explain the allowable service disciplines in detail 

• First Come First Served (FCFS) Under FCFS scheduling customcis aic 
served in the order in which they arrive The allowable service tinii 
distribution is negative exponential The service rate can be state dependent in 
which case )i,(j) will be the service rate at queue / when there are j customers 
in that queue In modelling a computer system an example use of the FCFS 
queue would be to represent a disk containing user files for a number of 
classes Since the basic operations performed at the device by the vinous 
classes are the same it is reasonable to assume that the average service times 
across classes are nearly equal The actual number of file accesses for 


customers of each class can be represented in the model by appropriate values 
of the visit ratios of each class at each queue 

• Last Come First Served with Preemptive Resume (LCFS PR) Under this 
discipline an arriving customer preempts the customer in service (if any) and 
Its service is begun immediately When a customer completes service the 
most recently preempted customer resumes service at the point at which it wis 
interrupted The service time distributions should have a rational Laplace 
transform This kind of service discipline can be used to model a CPU in i 
system where high priority system tasks are dispatched more frequently 

• Processor Sharing (PS) This is the idealization of round robin (RRi 
scheduling Under this discipline the control of the processor circul Ucs 
among all jobs in the queue Each customer receives a quantum of service 
before it must relinquish control to the next customer in the queue ind rejoin 
the queue at its tail The service time distributions should have a ration il 
Laplace transform This kind of service discipline can be used to model a CPI 
scheduling in multiprogrammed computer systems where some form of round 
robin scheduling is employed 

• Infinite Server (IS) This is also called dela} service discipline There will he 
simultaneous service for all the customers in the queue b\ an indnidu il slivci 
and hence the total delay in such a queue is equal to the service time at the 
queue Again the service time distributions should have a rational Laplace 
transform 

It has been explained earlier that FCFS queues are allowed to have a state dependent 
service rate An IS queue can also be thought of as a queue with state dependent service 
For example in a single class queueing network with a service rate of |i for each 
customer the IS queue will have a service rate of n\i when there are n customers m it 
The most important use of a queue with state dependent service is in the representation of 
Flow Equivalent Service Centres (FESCs) that can be used to reduce certain p irts of llic 
network and represent them as a single queue in a smaller possibly simplci queueing 


12 



network representation of the system The FESCs can also be studied in isolation In the 
next section we will discuss the compact representation of portions of a queueing network 
using FESCs 

2 3 Compact Representation of Portions of a Queueing 
Network 

Divide and Conquer is a well known technique to solve complex problems In electric il 
circuits by applying Norton s theorem any passive electrical circuit between two 
terminals can be replaced by an equivalent current source and a resistance in parallel to 
the current source The value of the current source is equal to the current that would flow 
when the two terminals are shorted in the original network The parallel resistance is 
determined by determining the resistance of the network when looking into it from the 
terminals This is illustrated m the example in Figure 2 3 We can study the Llectiical 
network shown in Figure 2 3a for various values of R by replacing the rcmiiniiu 
components of the circuit by a single current source whose current capicilv 1 is 
determined by shorting the points A and B (desired resistor R) and finding the curreni 
through this short circuit This is shown in Figure 2 3b The parallel resistance Rp is 
determined from Figure 2 3c Thus the electrical network of Figure 2 3a can be replaced 
b> the equivalent shown in Figure 2 3d 


R 

'b 

Figure 2 3 Norton s reduction in electrical circuits 


© 




-'WV“ 


a Original Network 




n 




d Equivalent Network 

Figure 2 3 Norton s reduction in electrical circuits (contd) 









IChandy Herzog and Woo described a method similar to Norton s theorem for 
reducing queueing networks They describe a technique to reduce a queueing network 
with N FCFS exponential service queues and M customers in the network In this 
technique a smaller network in which all queues except those in the subnetwork arc 
replaced by a single FES They show that for certain classes of system parameters the 
behaviour of the equivalent network is the same as that of the original network This 
result is called the Chandy Herzog Woo Theorem The technique is illustrated below 

Consider the queue Ql between points A and B of Figure 2 4a Obtaining the queueing 
statistics at this queue for different characterisations of it will be simplified if the rest ol 
the network can ha\e a more compact representation preferably as a single queue To 
achieve this we can construct an equivalent queue representation ot the rest of the 
network The method is known as Norton reduction and is described below The 
subnetwork is replaced by a state dependent queue with service rate p(j) which is i 
function of the number of jobs j in the queue This service rate p(j) is the same as the 
throughput T(j) between the points A and B of Figure 2 4b when there are j customeis in 
the network Note that the network of Figure 2 4b is the same as that in Figure 2 4a with i 
short across the queue Ql The equivalent network representation is shown in Figure 
2 4c In general the ponion of the network that can be isolated is not limited to a sinek 
queue and can be a network of queues Such a network of queues is called the designated 
netwrok and the remaining network is called the aggregate The aggregate is reduced to i 
FESC and the performance of the designated network is studied for the variation of the 
parameters of the designated network 

The hierarchical decomposition described above produces exact results for the large class 
of BCMP networks The decomposition method based on the Chandv Herzog Woo 
theorem is summarized below 

‘K M Chand} U Herzof^ and L Woo Parametric Aimhus of Qiteueing Networks IBM Jour ofRi-s 
and Develop 19 36 Jan 1975 

k M Chandj U HtrzoaandL Woo Approximate Anahsts of General Qiieueim, Vetwarki IBM 
Jour of Res and Develop 19 36 Jan 1975 
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• Select the designated subnetwork of queues from the original network that 
needs to be studied in detail The remaining portion of the network is the 
aggregate and will be reduced to a single RESC 

• Generate a queueing network in which the service times at all the queues m 
the designated network is set to zero i e they are shorted Note that the 
designated network should be selected such that the throughput through all the 
shorts in the new network obtained above should be identical 

• Solve the above queueing network using known techniques Solve this for all 
possible values of the network population in the original network The 
throughput through the short for different populations corresponds to the 
service rate of the FESC with that many in its queue 

• The FESC for the aggregate is now available 

• Obtain the equivalent network in which the designated network is m pi ill 
and the aggregate is replaced by its FTESC The queueing statistics ot the 
designated network in the equivalent network will be identical to those in the 
original network 

Aggregation produces exact results for queueing networks with a product lorm solution 
However aggregation of product form networks mav turn out to be computationally more 
expensive than simple MVA technique if the purpose is to solve tor onK one set ot 
parameters of the designated network The real advantage of aggregation is in solving non 
product form networks The components that do not have product form are put in the 
designated network and a FESC is obtained for the aggregate which contains onlv 
component queues that yield a product form solution The equivalent can be solved bv 
using approximate non product form solution techniques or by simulaiion Since solving 
an actual non product form network is computationally very expensive the reduced 
number of queues in the equivalent model would reduce the computing time However m 
this case assegation is only an approximation This is because the behav lour of the 
aggregate cannot be exactly replaced by the FESC since the information regarding the 
location of the customers in the aggregate is discarded However in many cases ot 
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practical interest it turns out that the results of the actual network and the equivalent 
network are close enough and the aggregation is a useful approximation 


2 4 Networks of GI/G/m Queues 

So far we have been dealing with queues with Poisson external amvals and service time 
distributions that were either exponential or had a rational Laplace transform Further it 
was only the mean service time that was used m the analysis We now consider queues 
where the external arrivals do not necessarily form a Poisson process and can ha\e i 
general distribution for the interarrival times We also allow the service times to be 
characterised by a general distribution An example of a network of such queues would be 
a packet switched communication network where packet arrivals form a burst> piocess 
not necessarily Poisson and the service times for each packet are fixed Such netW'Oiks 
can be modeled as a networks of GI/G/m queues The theory and alj^ori thins loi open 
networks of GI/G/m queues has been developed by Whitt^' ' This algorithm requiics th it 
the first two moments of the interarrival time and service time distributions bt known 
This technique can be used to solve an open network of queues having FCFS servicL 
discipline and multiple servers for single and multiple classes of customers Thu 
fundamental assumption made in this algorithm is that queues in the network are 
stochastically independent and a product form solution can be given The routing m itnv 
is used to obtain the two moments of the service time and interamval time distributions 
into the individual queues Each of the queues is then treated as an independent Gl/G/iii 
queue and a product form solution for the network is obtained 

A customer that is in the network can experience any of the following three b isie 
operations merging or superposition thinning or splitting and the departure of the 
customer from the network A multiplication factor associated with a queue y deteimines 


'W Whitt The Queueing NenwrkAnahser BSTJ Vol 62 No 9 Nov 1983 
W Whitt Performance of ihe Queueing NetiyorkAnahier BSTJ Vol 62 No 9 Nov 1983 
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if the customers are to be split (y > I 0) merged (0 0 < y < I 0) or not modified (y = 10) 
at that queue The routing matrix determines the next destination of a packet after its 
service at a queue 

In the following we briefly summarize the theory of the technique used in the Whitt 
algorithm used to analyse a network of GI/G/m queues Variability of a distribution is 
defined as the ratio of the variance to the square of the mean for that distnbution The 
arrival process of customers to any queue is a point process In this technique wl will 
approximate the point process corresponding to the arrivals by an equivalent renew il 
process using two parameters rate (mean number of arrivals per unit time) and thu 
variability (defined above) The approximation for a stochastic point process by a renew il 
process is described by Whitt This approximation involves the decomposition of the 
mode! by replacing the component flows (point process) by an independent lenew il 
process The distribution of the renewal process is obtained by fittin^ a distribution to i 
few moments of the interval between two point arrivals Thus an\ point process is 
approximated by a renewal process characterized by the mean and squared coefficient ot 
variance of the distribution that is fit to the moments evaluated from the properties ol the 
point process (Squared coefficient of variance is the variability parameter ) 

The Whitt algorithm converts the resulting piocess due to mer^in^ splittiiu and 
departure of traffic at the queues to a suitable renewal process (nett amv il pioeess) 
characterized by its mean and squared coefficient of variance This conversion is done b\ 
solving a set of linear equations for each parameter The queueing network is then 
decomposed into individual queues which are analysed separatelv Each queue is 
characterized by a standard GI/G/m queue with four parameters the mean and the 
squared coefficient of variance of the arrival process and the service time distributions 
Approximate two moment formulae are available to compute the performance of these 
queues in literature The def'omposition into independent queues can be interpreted as i 


'W Whitt Approximating a Point Process b\ a Renewal Procea I Two Basic Mitlwds Opurilmns 
Researeh Vol 30 No 1 pp 125 147 Jan Feb 1982 



generalization of the product form solution for the Markovian networks In Markovian 
models the population vector at each node are stochastically independent so that the 
probability mass function of the vector is a product of the mass functions of tht 
individual components of the vector The Whitt algorithm is desenbed in more det ul in 
the next chapter 

2 5 Fork- Join Queues 

There are situations when a given job is decomposed and distributed to be serviced in 
parallel and are aggregated after the sersice completion of the decomposed jobs Sucli 
queues are called fork join queues where the fork process corresponds to the 
decomposition and service in parallel and the join process corresponds to the aggregation 
after completion of all the sub jobs A fork join queueing system contains multiple 
parallel queues and can operate wither without synchronizing queues A synchronism^ 
queue is a waiting place for each of the parallel queues (sibling queues) where the 
customers wait for aggregation after receiving service at that queue It theic is no 
synchronising queue the customers will wait at the server and block the server tor other 
customers in that sibling queue In general the sibling and s>nchromsing queues in i 
fork join node will have infinite queueing capacities use FCFS service discipline anil 
operate independentlv fork join queue is illustrated in Figure 2 5 ^n incoming job is 
split at forking node A and aggregated at the joining node B Even, job arriving it the 
forking node consists of the same number of subtasks as the number of sibling queues 
and upon arrival at the forking node is immediately divided into a number of subtasks 
and each is simultaneously submitted to the sibling queues 

An exact analysis of a fork join queue is difficult However many approximate models 
have been described in literature and we have used the techniques described in [’ '] 

'C Kim and A K Agrawala Aiuih sis of the Fork Join Queue IEEE Trans on Compuicrs Vnl 18 

No 2 pp 250 255 Feb 1989 

Y C Liu and H G Perros A Decomposition for the Anahsis of a Closed Fork/Join Qneiitini, 

S\stem IEEE Trans on Computers Vol 40 No 3 pp 365 370 Mar 1991 
’Y C Liu and H G Perros Approximate Ancihsis of a Closed Fork/Join nwdet European Journ il of 

Operlational Research 53 pp 382 392 1991 




A fork join queue may or may not have a synchronising queue If there is no 
synchronising queue the subtasks wait at the server and they are rejoined after the I isi 
task has finished its service from the sibling queues Thus if there are k sibling qiiciits 
they can be replaced by an equivalent queue whose service time distribution is equal to 
the distribution of the maximum of the k service time distributions of the sibling queues 
For QNAT we have assumed exponential service time distribution at the siblings noi 
necessarily identical If the network is open we replace the fork join queue by in 
equivalent queue whose service time distribution is the distribution maximum of the k 
exponentials K the network is closed we first reduce the fork join queue to a single 
equivalent FCFS queue with the service time distributior given by the distribution of the 
maximum of k sibling service times We then fit an exponential distribution to this 
distribution that has the minimum mean square error Fork join queues without 
synchronising queues are illustrated in Figure 2 6 

A fork join queue usually has a synchronising queue so that the subtasks can le ive the 
sibling servers asynchronously The subtasks after completion of service at the sibling 
queues wait in the synchronizing queues When all the subtasks finish service theN 
recombine and leave the queue In this kind of a queue the time spent in the 
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a A fork join queue without synchronizing queues 



b Equivalent reduced queue 

I 

Figure 2 6 Fork Join without synchronizing queues 

synchronising queue waiting for completion of the subtasks can be a significant portion of 
the response time (The response time for a job is the time between the forking and the 
joining instants of the job ) Liu and Ferros have developed an algorithm for the analysis 
of a closed queueing network with synchronizing queues The algorithm is an iterative 
approximation and for a k sibling fork join queue with M customers it works as 
follows Obtain a two sibling closed fork join queue by choosing any two ot the siblings 
from the original k This fork join queue with j customers can be solved to obtain the 
throughput for j=\ M This is used to obtain the Norton equivalent for the two siblin_ 
fork join queue Another two sibling queue is now obtained in which one of the siblings 
is the Norton equivalent from above and the other is one of the remaining siblings in the 
ongmal k sibling queue This can be solved as before and the process repeated till all the 
it siblings have been used to obtain the Norton equivalent of the original k sibling fork 
join queue A more detailed description of the algorithm follows and the implementation 
details are given in the next chapter 





Synchronizing 

queues 


a The fork join queue 



b Equivalent reduced queue 
Figure 2 7 Fork Join with svnchronozing queues 


Consider the closed queueing network with M customer and a k sibling fork join queue 
shown in Figure 2 8a To apply Norton s theorem between points A and B of Figure 
2 8a the network shown in Figure 2 8b is used To apply Norton s theorem consider 
sibling queues 1 and 2 as shown in Figure 2 8c This system can be sohed exacth as 
follows Let P,j be the probability that there are i jobs in sibling queue 1 and j jobs in 
sibling queue 2 This is a closed queueing network with state dependent sersice rates at 
the queues The state transition diagram forM=3 is shown in Fig 2 9 In this figure (/j) 
represents the state of the two queue system and corresponds to having i customers m 
sibling queue 1 and j customers in sibling queue 2 It is easy to see that states (0 2) (1 2) 
(2 2) (0 1) (1 1) (2 1) (0 0) (10) and (2 0) are not possible Thus the total number of 
feasible states are ( V/+ 1 )^ M’ 

From the state transition diagram of Figure 2 9 for A/ = 3 we have 
Pv/ \t ( )ii(Af) + |J.->(A0 ) = Paiw i P-i(^ + Pw ijW li-’(Af) 


2 ^ 



( ^ll(A/) + H^Cn) ) = 1 Hi(A/) + PMn^\ (^+0 

Pmq^i(A0 = Pmi 1^2(1) 

^nM ( + M-t ( M) ) = Pn+\M Hl(«+1) + P/I IjM 

Po^ \y-2iM) = Pi^ |ii(i) 

M-\ U-\ 

Pa/^+ X Pa/.+ 2; Pj>/=10 

=0 ;=o 


Solving this set of equations we get 

w 

P/ijW = Pww J~J P'’(7) Q<n<M \ 

y= I 

M 

Pwi = Pa/a/ Pi(7) 0<n<M \ 

j= i 


Pw vy = 


w-i 


' - I 

sO 

where 


1 0 



P "(;) + 


M I 

s 


i 0 


M 

n 

/=i. + I 


Pi(0 


Pi (y) = 


^ l (;) 

Pi(M) 


p (;) = 


Pi (; ) 
p 1 M ) 


The system throughput of Figure 2 8c is given by 

M 1 w-i 

Ueff(M)= pi(M) X P>iV^ + Pi(^^) X 

=0 y-O 


l<n<M I 
(2 8 ) 
1 <n<M 

(2 9) 


(2 10 ) 

(2 11 ) 


(2 12 ) 


24 



b The k sibling fork join queue as a closed queueing network 

li|{ i) 



c Siblings queues 1 and 2 are reduced to an FESC with rate |(j) 


Figure 2 8 Fork join with synchronizing queues for closed network 




0 |( 1 ) 



1=1 2 
j=12 

a + b = M 


d Sibling queue 3 and FESC from the reduction of siblings 1 and 2 are 
used to obtain the FESC for a three sibling fork join queue 
The service rate of this FESC is •> 

Figure 2 8 Fork join with synchronizing queues for closed network (contd) 



This Norton equivalent can be used with a third sibling queue of the k sibling queue to 
obtain another equivalent queue and this is shown in Figure 2 8d Thus by applying this 
technique iteratively k 1 times the k sibling queue finally reduced to a single queue with 
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state dependent service rate The implementation details of this algorithms is given in the 
next chapter 


2 6 Conclusion 

In this chapter we have explained the theory and the outline of the various algorithms th it 
have been used in QNAT We started with a brief introduction to product form solutions 
to open and closed queueing networks Calculation of the normalising constant is the 
first bottleneck in the use of the closed queueing networks in modeling practical systems 
A convolution algorithm to calculate the normalising constant was gi\en Mean value 
analysis a technique developed to overcome the numerical difficulties m the calculation 
of the normalising constant was then described Multiple classes of customers can be 
handled through BCMP networks We also described techniques to reduce parts ot i 
large complex queueing and represent it by a flow equivalent server Fork join queues 
and approximate algorithms to anal>se them were also described We also discussed 
networks of non Markovian queues 

In the next chapter we will outline the algorithms that we have used in QNAT to solve 
BCMP queueing networks networks with fork join nodes and open netwoiks of non 
Markovian queues We also give examples to use QNAT and explain the input 
mechanism to QNAT 
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Chapter 3 


Methods Used for Solving Networks with Infinite 
Capacity Queues 


In this chapter we describe the algorithms used in the QNAT package and provide 
examples to illustrate the use of the package and also to venfy the working of the 
algorithm implementation For each type of analysis technique available through QNAT 
we will give a brief introduction to the technique descnbe the input and output 
parameters present the algorithm and then show some examples to explain the use of the 
method and also to verify the correctness of the algorithm implementation 


3 1 Queueing Networks with a Single Class of Jobs 

In these cases there will only be one class of customers (jobs) in the network i e all the 
customers in the network are treated to be of the same class The network can be either 
open or closed for this class of customers 

311 Open Net>vorks 

In open networks the jobs come from the externa! world (i e outside the network) into 
the network and after getting service at different queues i e nodes in the network the\ 
depart to the external world In an open network there will be at least one node where 
jobs from the external world will enter the network The algonthm used for the analysis 
of an open queueing network with single class input is the GI/G/m method ' This 
method requires the first and second moments of the interamsal time and service time 
distnbutions It computes the first and second moments of the net amval process at each 


'w Whitt The Queueing Network Anahser BSTJ Vol 62 No 9 No\ 1981 
W Whitt Performance of the Queueing Network Aiuihser BSTJ Vol 62 No 9 Nov 1981 
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node in the network by solving a set of linear equations The performance of the network 
IS then obtained by analyzing the individual queues in isolation Each queue is a standard 
GI/G/m queue characterised by the first two moments of the arrival process and the 
service times Details on this method may be obtained from 

Assumptions 

• All the queues in the network have infinite buffer capacity 

• The service discipline at all the queue in the network is FCFS 

• There can be creation or combination of jobs at the queues m the network 

• The routing matrix for the jobs in the network is static 

Input Parameters 

N number of nodes in the network 
m, number of servers at node z for i = 1 2 M 

Xoi mean of the external arrival rate distribution at node z for z = 1 2 N 
Coi^ squared coefficient of variance of the interarrnal time distnbuiion at node z 
for z = I 2 aV 

X, mean of the service time distribution at node z for z = 1 2 -V 
Cu squared coefficient of variance of service time distnbution at node z 
for z = 1 2 M 

y, multiplication factor at node z for z = I 2 
[P] p,j routing matrix for z y = I 2 N 

variance 

C Squared coefficient of variance (SqCVar) ^ 

mean" 

The multiplication factor y, at node z specifies the creation or combination of jobs at that 
node If yi > 1 then jobs are created if y, < 1 then jobs are combined and if yi = 1 then 
jobs are neither created nor combined 


' M N Umeih A VV;«^/ou based Package for Queueing Nen^o^k Anahsis M Tech ihcMs IIT Kanpur 
Feb 1997 
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If Xj IS the nett arrival rate at node j then departure out of node j will be Xj Yj For 
example consider a switch in a packet communication network After processing a 
packet at the switch the packet can be split into a number of smaller packets and then 
transmitted to other nodes in the network At the destination these packets can be 
recombined to get the original packet 

Output Parameters 

X, mean nett arrival rate at node i for i = 1 2 ^ 

Cai squared coefficient of variance of nett arrival rate distribution at node i 
for I = 1 2 M 

EN, mean number at node i for i = 1 2 

CnI' squared coefficient of variance of number at node z for / = 1 2 jV 
EVV, mean waiting time at node i for / = 1 2 

C, , squared coefficient of variance of waiting time at node / for z = 1 2 N 
dep, departure rate at node z for z = 1 2 JV 
V, visit counts at node z for z = 1 2 JV 
Ui utilisation at node z for z = 1 2 N 


Algorithm 3 1 GI/G/m solution technique for open networks ’ 


T, 

Ti - 

1 Pn 


for z = 1 2 

N 

^S1 ~ C'si + Pii 

(1 Cs.’) 

for z = 1 2 

N 

II 

r 

+ 

1M< 


for z = 1 2 

N 

qij - Yi Pij “ 

amval rate from node i to node; forz;=12 N 


'w Whitt The Queueing Net\wrkAnal\ser B S T J Vol 62 No 9 Nov 1983 
M N Umesh A Window based Package for Queueing Network Anahsis M Tech thcMs IIT Kanpur 

Feb 1997 
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( Algorithm 3 1 continued ) 


where 


X,x, 


[/.= 


for / = 1 2 N 


m. 


Cai = |ii + X ^ 0.1 Vj, for I = 1 2 N 


j 1 


^1, = 1 + CDi [{ goi Co,’’ 1)+X Pji)+^jPj. c/xj}] for/=12 N 


j I 


V)! ~ ^ Pji 9ji 1 Cj ) for i J ~ \ 1 N 


tj = 1 + /?ij ” [ max (Csj 0 2)1] for; = 1 2 N 

1 0 


rui = 


(1+4(1 D) 


for ! = 1 2 N 


1 0 


\)j=‘ 


for; = 1 2 N 


r-l 


qji 


<?j.= 


i V 

i-l 


for / ; = 1 2 N 




£N. = C, /?i, + £W, 


for f = 1 2 N 


M + C, C Pq^ \ 

C, „' = /»,[/. + f/, PJ +'«. fori = 12 N 

(I C.)' ^ 
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( Algorithm 3 1 continued ) 

'C. V, {Cj + ) p, 

y;/ oi/G/i _ 

2(1 V,) 


for I = 1 2 


2(1 f/,)(l X')’ 

( ) ^ j 

e 3 t/, (Ca,' + C,’ ) 

* ^ otherwise 


C., =’ 


Pq, (2 Pq, )T.' 

[//!,( 1 (/,)]' By, 


for 1 = 1 2 N 


for j = 1 2 


n 

dep, = X, y, { 1 X Pu ) 


for 1 = 1 2 N 


for 1 = 1 2 N 


2 ^)j 


Network Performance 


II 

Throughput = J 

/-I 

V 

Total departure rate = ^ dep, 

1=1 

V 

Mean number = ^ £A( 


Squared coefficient of variance(mean number) = 


S ( C,. £N, ) 


(Z 


n 

Sojourn time = Z ^ 
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Example 3 1 


X=0 45 

d 

^p=0 5 

— ► 

"'S 

p=0 5 

— ►(!) 


^(2) 



p=l 0 


Figure 3 1 The open network for Example 3 1 


Description of the network for Example 3 1 

We have a two node single class open network There are external arrnals onl\ to node 1 
The interarrival time distribution for this is hyperexponential with mean 0 45 and SqCVar 
2 25 Node 1 has h\perexponential ser\ice time distribution with mean 1 0 and SqCVar 
2 25 Node 2 has Erlangian 4 service time distribution with mean 1 0 and SqC\ ar 0 25 

Input parameters to Example 3 1 



Node 1 

Node 2 

;n, 

1 

1 


0 45 

00 

c,,' 

2 25 

1 0 

T, 

1 0 

1 0 

c.' 

2 25 

0 25 

fl 

1 0 

1 0 


Table 3 1 Input parameters of Example 3 1 



Node 1 

Node 2 

Node 1 

00 

05 

Node 2 

1 0 

00 


Table 3 2 Transition probability matrix of Example 3 1 
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Output pirameters forExamnle 1 from QNAT 



Node 1 

Node 2 

^)l 

0 45 

00 

Co,’ 

2 25 

1 0 


0 90000 

0 45 

Ca.' 

1 749342 

1 577437 

£N, 

17 09734 

0 786415 

c,: 

0 3078830 

2 405392 

E\V, 

17 99704 

0 747588 

c 

0 3056561 

4 12567 


0 45 

0 

V, 

2 0000 

1 0000 

c, 

0 9000 

0 4500 


Table 3 3 Output parameters of Example 3 1 


3 12 Fork Join Analysis (Open Networks ) 

The performance of distributed and parallel processing systems can be modelled by fork 
join analysis In this analysis the sibling queues can be with or without synchronizing 
queues If there are k sibling queues then each job will be split to k subjobs and thi. 
subjobs are served independentlv of one another The jobs will recombine when all ths. 
subjobs complete service in the sibling queues In the case of fork join queues without 
svnchromzing queue the jobs will recombine when the last subjob finishes its service at 
Its sibling queue In other words recombination of the subjobs will be followed by the 
departure of the subjob coming out of the sibling queue that has minimum service rate 
1 e the subjob spends maximum time in that sibling queue Thus the effective service 
time distribution of the sibling process will be maximum of the k random serv ice times 


Assumptions 

• All the sibling queues will have infinite buffer capacity 

• The service discipline at all the sibling queues are FCFS 



















• There will be no external arrivals into the sibling queues (Xo,) 

Xo, = 0 0 for / = 1 2 k 

There will be neither creation nor combination of jobs at the sibling queues 
The multiplication factor (y) is 1 0 at the resultant single queue 7=10 

• All the sibling queues are single server queues 

• The service time distribution at all the sibling queues are exponential 

Input parameters 

k number of siblings in the fork-join queue 

l/|li mean service time at sibling queue t for i = 1 2 L 

Output paiameters 

T mean resultant service time distribution of the fork join model 
Cs squared coefficient of variance of resultant service time distribution 

Algorithm 3 2 Fork join without synchronizing queues for open networks 


1 /^ 1 , fori = 1 2 k 

p, X 

F,, ( X ) = 1 e (Exponential distribution with mean |i, tor queue i) 

for / = 1 2 k 

Pi X 

/si ( X ) = Pi e (Density corresponding to above distribution) 

for z = 1 2 k 

n 

Fw(x) = n F 1 ( X ) (Service time distribution = maximum of k exponentials ) 
1 

K K 

/w(x)=X (/si(x) n Fsj(x )) (Density corresponding to above distribution ) 
(=1 } 1 
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( Algorithm 3 2 continued ) 


= J /w ( X ) X dx (Mean of the Fw(x) ) 

0 

'I" = I /w ( X ) x^ dx (Second moment Fw(x) ) 

0 

= j (Squared coefficient of variance of Fw(x) ) 

(T)’ 



Description of the network for Example 3 2_ 

There are two sibling queues without s>nchronizing queues in the sibling process The 
mean service rate at sibling queue 1 is [i| and in sibling queue 2 is p. as shown in Figure 
3 2 At both the sibling queues the service time distribution is exponential The 
equiv alent reduced network is then obtained as shown in Figure3 2b This open network 
can then be solved using Algorithm 3 1 
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Input parameters to Example 3 2 



Sibling 
queue 1 

Sibling 
queue 2 

1/p, 

20 

10 


Table 3 4 Input parameters of Example 3 2 


Output parameters for Example 3 2 from ONAT 



Resultant 

queue 

T 

23 33333 

c; 

0 673469 


Table 3 5 Output parameters of Example 3 2 

For a fork join queue with synchronizing queues after completion ot service at the 
sibling queues the jobs will stay in the synchronizing queues When all the subjobs finish 
sen ice only then are the jobs recombined Analysis techniques for such queues are not 
asailable in the literature Therefore this model has not been implemented in QNAT 

3 13 Closed Networks 

A t-losed network is one in which jobs are circulated inside the network There will be no 
external arrivals to any of the nodes or external departures from anv of the nodes in the 
nenvork The algorithm used for the analysis of such closed networks with single class 
input is state dependent MVA'‘ This is an iterative algorithm which requires only tht. 
mean service time (first moments) of the service time distributions 

Assumptions 

• All the queues in the network have infinite buffer capacity 

• There are no external arrivals to any of the nodes in the network 


* Edward D Lazowska JohnZahorjan G Scott Graham and Kenneth C Sescik Quantitcitne S\snni 
Performance Prentice Hall Inc 1984 

Michael K Molloy Fundamental of Performance Modelling Macmillan Publtshma Company 1989 





• The routing probabilities for the jobs in the network are static 

• The service time distribution of the queues in the network are exponential 

Input parameters 

N number of nodes in the network 

m, number of servers at queue i for j = 1 2 ^ 

T, mean service time at queue i for i = ] 2 N 
QUpe^ service discipline at queue i fori=12 ^ 

FCFS LCFS PR PS and IS 
[P] p,j routing matrix 17=12 ^ 

Output parameters 

mean number at queue 1 for i = 1 2 N 
EWi mean waiting time at queue i for / = 1 2 N 
A, mean throughput at queue i for i = 1 2 N 
Vi visit ratio at queue i for 1 = 1 2 N 


Algorithm 3 3 State dependent MV \ solution technique for closed networks 


M 

number ot customers m the closed network 



N 

number of nodes in the network 




tij 

mean serv ice at queue i when j in the queue for i = 

1 2 

N 




for; = 

1 2 

M 


= t, if Qr\pe i IS IS 















if J < m, 






f 

J 

for I = 1 2 

N 





for; = 1 2 

M 




x, 


Qt\pe i = FCFS LCFS PR or PS 




if J > m, 






V tn, 






39 



( Algorithm 3 3 continued ) 


The visit ratios are obtained by solving the set of linear equations 
V, = X Pji for I = 1 2 


j=i 

Initialisation 

N,(0) = 0 for r = 1 2 ^ 

jmp,(0 0) = 1 for I = 1 2 N 

jmp,0 0) = 0 for 7 = 1 2 M for z = 1 2 ^ 


Do steps 1 2 and 3 for c = 1 2 
Step 1 


M 


w 


.(c) = X j''^.(j)jmp,(j 1 c 1) forz=12 

N 

W(c) = X ^1 w,(c) 


N 


X{c) = 


W{c) 


Step 2 

N,(c) = Xic) V, w,(c) for / = 1 2 
Step 3 

jmp,(j c) = >dc;) V, T| jmp,0 I c) for z = 1 2 jV and for; = 1 2 


jmp,(0c)=l X jmp.O c) for z = 1 2 




;-i 


£jV, = N,(Af) By, = w,(A/) X, = X(M) V, for z = 1 2 N 



40 





Description of the network for Example 3 3 

There are 3 nodes and 6 jobs in the network Node 1 uses a PS service discipline Node 2 
and Node 3 use the FCFS service discipline This example is taken from Molloy*'' This 
IS a model of a simple computer system with Node 1 as the processor and Nodes 2 and 3 
as the I/O devices The service time distributions at all nodes are exponential The mean 
service time at Nodes 1 2 and 3 are 0 028 0 04 and 0 28 respectively 


Input parameters to Example 3 3 



Node 1 

Node 2 

Node 3 

fn, 

1 

1 

1 

X. 

0 028 

0 040 

0 280 

Qr\pe, 

PS 

FCFS 

FCFS 


Table 3 6 Input parameters of Example 3 3 



Node 1 

Node 2 

Node ^ 

Node 1 

0 1 

07 

02 

Node 2 

1 0 

00 

00 

Node 3 

1 0 


00 


Table 3 7 Routing matrix for Example 3 3 


Output parameters for Example 3 3 from QN^ 



Node 1 

Node 2 

Node 

E\ 

0 886640 

0 886640 

4 22672 

BV, 

0 02310 

0 03300 

0 93800 

h 

17 35107 

12 14575 n 

3 47021 

V, 

1 00000 

0 70000 

0 20000 


Table 3 8 Output parameters of Example 3 3 


M,.haei K Mollm Funjamma! of Modelling Macmillan Publi.hm Companj NS!) 
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314 Fork Join Analysis (Closed Networks) 

ase of a fork join queue without synchronizing queue in a closed queueing 
network the jobs will be combined after the last subjob has finished the service at its 
sibling queue In this the reduction step is a two fold one The first step is the same as the 
one stated in Section 3 1 2 for the case without synchronizing queue i e getting the 
resultant distribution of the service time In the second step we fit an exponential 
distribution to this resultant distribution such that the mean square error between the two 
distributions is minimum 

Assumptions 

• All the sibling queue have infinite buffer capacity 

• The ser\ ice discipline is FCFS at all the sibling queues 

• The ser\ ice time distribution at each sibling queue is exponential 

• Each sibling queue is a single server queue 

• There are no external arrival to any of the sibling queues in the model 
Input parameters 

k number of sibling queue in the fork join model 

l/(i, mean sen ice time at the sibling queues / = 1 2 k 


Output parameters 

T mean resultant ser\ ice time of the fork join model ( exponential fit ) 


Algorithm 3 4 Fork join without synchronizing queues for closed networks 


l/|li for 1=12 

k 


Fs, ( X ) = 1 e 

(Exponential distribution with mean p, for queue i) 
for 1 = 1 2 k 

p, X 

Ai ( X ) — Pi G 

(Density corresponding to above distnbution) 
for 1 = 1 2 k 
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( Algorithm 3 4 contmued ) 

K 

^w(x)= n ( X ) (Service time distribution = maximum of k exponentials) 
/■=! 

Clx 

f cstm ( X ) = 1 e (An exponential fit to Fw(x) ) 

M 

error (X C)= | [Fw(x) F„im(x)]’dx 
0 

Omin = mm [ error ( x U) ] 

Cj 


t — 1/ Umin 


Example 3 4 




b The equivalent reduced network 

Figure 3 4 Fork join queues without synchronizing queues in a closed network 
for Example 3 4 ____ 
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PescnEtio n of the network for F.^amp lp ^ a 

There are two sibling queues in the sibling process There are no synchronizing queue for 
either of the sibling queues The mean service rate at sibling queue 1 is n, and sibling 
queue 2 is p. At both the sibling queues the service time distribution is exponential 
The resultant service time distribution is also exponential The equnalent network m 
Figure 3 4b IS a closed network and can be solved by using Algorithm 3 3 

Input parameters to Example 3 4 



Sibling 
queue 1 

Sibling 
queue 2 

1/p, 

20 

10 


Table 3 9 Input parameters of Example 3 4 
Output parameters for Example 3 4 from ONAT 



Resultant 

queue 

T 

23 64108 


Table 3 10 Output parameters of Example 3 4 

In the case of fork join with synchronizing queues after completion of service at thi. 
sibling queue the subjobs will wait in the synchronizing queues When all the subjobs 
finish their service the\ are recombined The solution procedure for such a queue is 
based on decomposition and aggregation This algorithm is from ' There are two 
procedures to analyse this situation One is a numerical procedure which can be applied for small 
s\ stem ( = 2 sibling queues ) In this procedure the following steps are carried out 

• Generation of all states of the system and the corresponding rate matrix Q 

• Solution of the linear system iQx^ = 0 steady state probability distribution 


' "i C Liu and H G Perros A Decomposition for the Anahsis of a Closed ForUJoin Qiteming 
System IEEE Trans on Computers Vol 40 No 3 pp 365 370 Mar 1991 

\ C Liu and H G Perros Approximate Anahsis of a Closed Fork/Jom model European Journal of 
Operational Research 53 pp 382 392 1991 
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The second approach is an approximate procedure which can be applied for large systems 
2 sibling queues) For such networks we first obtain the FESC for the fork join queue 
by shorting the input and the output of the fork join node This FESC is then used in the 
larger network which is solved using techniques described before The FESC for the fork 
join queue is then obtained as follows 

• Two siblings are selected and solved exactly to obtain their FESC 

• This FESC is used to successively aggregate the other k 2 siblings into the 
FESC 

Assumptions 

• All the sibling queue have infinite buffer capacity 

• The service discipline is FCFS at each of the sibling queues 

• The service time distribution at each sibling queue is exponential 

• Each sibling queue is a single server queue 

• There are no external arrival to any of the sibling queues in the model 

Input parameters 

k number of sibling queues in the fork join model 
M number of customersQobs) in the network 

T,( j ) = l/pX j ) mean service time at sibling queue i whenj customers in it 

Output parameters 

"Ceff ( J ) resultant mean service time of the fork join mode with j jobs m it 

Algorithm 3 5 Fork join with synchronizing queues for closed networks 

Do steps 1 & 2 for c = 1 2 k ] 

Stepl Do (A) (B) and (C) for; =12 M 

(A) n 0<n<;l 

k- +1 




( Algorithm 3 5 continued ) 

~ - - - 

(B) Pj„ = n p,( k ) 0<n<j \ 

*= +1 

1 0 

(C) 

' + X n pc+i(<?) +x n pc(n 

-0 ti=+l I 0 I I \ 

where 

P<-+i( O lOcC / ) 

p ( , ) = for j = 1 2 j 

Step2 

^ ;-l 

fk-i-i (J)= [K.iij) X P.J + M.(2) X Pj 

-0 =0 


t.ff(;) = l/[pAi(7)]=l/[Ue,K;)] ;=12 M 


The GUI in QNAT allows us to input the service tune at the sibling queues as a table or 
as an expression For the case of table the mean service time at each sibling queue is 
entered separately for each value of the number of jobs at that sibling queue Example 3 5 
will illustrate a situation where the mean service time input is gi\en as a table For the 
case of an expression the mean service time is entered as a function of // the number of 
customers in the sibling queue The expression can be a constant or a function of the 
variable n in the standard mathematical form Example 3 6 will illustrate a situation 
where the mean serv ice time is input as an expression 

Description of the network for Example 3 5 

There are two sibling queues with synchronizing queues in the closed network with four 
customers in the network jobs at point A in Figure 3 5a will be split into two and passed 
on to the sibling queues After service at the sibling queues the subjobs will wait in the 
respective synchronizing queues When both the subjobs finish at their respective sibling 
queues they recombine at point B in Figure 3 5a The service rate at both the sibling 
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queues are state dependent and have an exponential distribution This service can be 
given in two forms either as a table or as an expression The equivalent network in 
Figure 3 5b is like a two node network with service centres having state dependent 
service rate This may be analyzed using Algorithm 3 3 to get the performance of the 
equivalent network 
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Example 3 5 Service time as a Tahl^ 
Input parameters to Example 


Number 

1 

2 

3 

4 


10 5 

10 25 

10 125 

10 0625 


5 66667 

5 44444 

5 296296 

5 19753 


Table 3 1 1 Input parameters of Example 3 5 


Output parameters for Example 3 5 from ONAT 


Number 

1 

2 

3 

4 


12 486 

10 683 

10 238 

10 092 


Table 3 12 Output parameters of Example 3 5 

Example 3 6 Service time as an Expression 
Input parameters to Example 3 6 



10 + 2 { n ) 


5+ 1 5'^( n ) 


Table 3 13 Input parameters of Example 3 6 


Output parameters for Example 3 6 from QNAT 


Number 

1 

2 

3 

4 


12 486 

10 683 

10 2^8 

10 093 


Table 3 14 Output parameters of Example 3 6 

3 2 Queueing Net\^orks with Multiple Classes of Jobs 

When there are more than one class of customers in the network the system can be open 
with respect to some classes of customers and closed with respect to the other classes If 
the number of closed class customers is zero than the system is an open network If the 
number of open classes is zero then the system is a closed network Otherwise the 
network is called a mixed network 


4S 
















In developing a queueing network model of a system we may need to allow a customer to 
change classes after service at a queue and be routed according to the routing probabilities 
of that class It has been shown that such a general model can be trivially mapped into a 
BCMP network that does not allow change of class Therefore in QNAT we do not 
allow customers to change classes 


3 21 Open Networks 

The network is considered to be an open network when all the classes m the network arc 
open In the case of multiple classses of jobs there must be external arrivals/departures 
to/from the network for each class This part has been earlier implemented by Umesh and 
the details of the algorithm and its implementation can be obtained from '' ' 


3 2 2 Closed Networks 

In a closed network all the jobs of all the different classes are to be circulated inside the 
system There will be no external arrnals to the system and hence no departures out ot 
the system The technique used to esaluate the closed queueing network is the Mean 
Value Analysis (MVA) technique The parameters like mean queue length mean waiting 
time mean throughput etc in closed multiple class queueing network which ha\e 
product form solution can be computed recursnely without computing the product terms 
and the normalizing constants Howe\er from this technique we can only get the first 
moments In spite of this limitation we have decided to use this algorithm because of the 
Simplicity of its implementation 


This algorithm is based on the Mean Value Theorem This theorem states that a 
customer of class k arriving to a queue in the network would see the same average 


'WardWhitt The Queiiemg NetywrkAnahser ® ^ 

M N Umesh A Windosy based Package for Queueing Nenvork Analysis M Tech thesis IIT 

' mTcZt S S Tuasenberg Mean Value Anahsis of Closed Multicham Queueing Neluork, 
ATM Vol 27 No 2 DP 313 322 April 1980 

■' Michael K Molloy Fundanienial of Performance Modelling Macmillan Publishing Company 


J of 
1989 
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number in the queue as an external observer would see if the network had one less 
customer of class k The concept of visit ratio is used in the algonthm Visit ratios are 
defined as the ratios of throughputs Visit ratio at node k Vk is the ratio of throughput of 
node k to the throughput of the reference node in the network For the analysis of closed 
networks node 1 is arbitrarily selected as the reference node in QNAT to calculate the 
visit ratios It is easy to see that Vk = 'ky/X[ where X, is the absolute arrival rate to node i 
In mean value analysis we first initialise and compute the waiting time at each node for 
each class in the network Next we compute the throughput of each class followed by the 
mean queue length at each node for each class in the network Finally we update the 
queue length distribution These parameters are computed iteratively until we reach the 
total population of each class 


Assumptions 

• All the queues in the network have infinite buffer capacits 

• There are no external arrivals to any of the queues for anv ot the classes 

• Classes are independent and the routing probabilities for eich class in the 
network are static 

• Exponential service times for all classes at all nodes in the network 

• Onlv the following four service discipline are considered in the network 
FCFS LCFS PR PS and IS and are indicated in the GLI as follows 

FCFS 1 LCFS PR 2 PS 3 IS 4 


• There is at least one class of jobs with a non zero population 
Input parameters 

N number of nodes in the network 

B total number of classes 

m, number of servers at node i for / = 1 2 

T,k mean ser\ ice time at node i of class Afori=12 // & k = I 2 J? 

Qtypei service discipline at node 1 for/=12 

C population vector ={C, C, C.) where a ts the population of class ^ 

[QJvvA qu.j) routing matrix for 1 7=12 ^&k=\2 R 
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Output parameters 


iSNik mean number of class k customers at node /forj=12 A^&^=12 R 
EWi k mean waiting time of class k customers at node i for i= 1 2 

k=l2 R 

X-,k mean throughput of class k customers at node j for r = 1 2 

k^\2 R 

V,k visit ratio of class k customers at node ffori=12 N & k = I 2 R 

R 

EN, = X EHv for 1 = 1 2 N 

t=i 


Algorithm 3 6 Multiple class MVA solution technique for closed networks 


n 1 


N number of nodes in the network 
R number of classes 

C {C| C Cy?) population vector where Ck IS total population of class A, 

0*=: { 0 0 0 ) zero vector zero in each class for all R classes 

^={00 10 0} unit vector in the k*'’ position 


1) Initialisation 

£^,( 0 ) = 1 0 for 1 = 1 2 

P,( 0 0 ) = 1 0 for 1 = 1 2 

PiiJ 2 ^ = 00 fori = 1 2 


N 

N 

N 


2) Do steps 3 4 & 5 for Cl = 0 1 2 C\ c->-0 12 

3) Do for A = 1 2 R &i=\ 2 N 

L 


c 


ck=0 12 Ck 


w,k( C ) = 


t,kV.kL/V(C) 

Ilk V,k 

1 

1 T,i. v,uiiv,(C)+ S ("'■ 

V 


Wi,( c ) = X "'a £ ) 


'M Re,seran<iS S La^enterg M™,, V,l„. Amhm oj Closed Qoeoeo,, •deou.rks 

ACM Vol 27 No 2 pp 313 322 April 1980 


> kTcLT moii‘;“ 

‘933 


1 





'Is®# 

2 % i' 

i ii V 


J of 
1989 
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( Algorithm 3 6 continued ) 

Clc 

MC)= 

Wk(C) 

ft 

4) LN,{C)=\ + Y, >.k(C) VaW,k(C) for/= 1 2 ^ 

k I 

5) Do for /= 1 2 JSI 8c j-M m, 1 


P.(; C)= ( I/j)* X ^k(C)V.kT.KP,(j 1 C-eu) 

i-i 

R -I 

X >^k(C) V„T., + X i) 

A-l J 1 

P,(0C)=1 

m, 


C {Q C Q } 

EN, = LN,{C) 1 for ;= I 2 .V 
EN,i. = k( C ) w,i,( C)for/=12 N Sc k=\ 2 R 
w,k ( C ) V,k T,k 

- forj=12 M Sc k - \ 2 R 

V,k 

?t.k(£) = MC) V.k for 1 2 N8ck=\2 R 


Example 3 7 

Description of the network for Example 7_ 

There are 2 classes and 3 nodes in the closed network There are a total of 6 customers 
with 3 in each class The service discipline for both the classes at node 1 is PS and at 
node 2 and node 3 are FCFS The mean service time for both the classes at node 1 is 
0 028 node 2 is 0 040 and node 3 is 0 280 This is similar to the one given in Example 

3 3 
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Figure 3 6 Multiple class Closed Network( 2 classes ) for Example 3 7 


Input parameters to Example 3 7 



Node 1 

Node 2 

Node 3 

7U| 

1 

1 

1 

Qnpe, 

3 

1 

1 


Table 3 15 Input parameters of Example 3 7 


T,k 

Node 1 

Node 2 

Node 3 

Class 1 

0 028 

0 040 

0 280 

Class 2 

0 028 

0 040 

0 280 


Table 3 16 Mean service time of Example 3 7 



Node 1 

Node 2 

Node 3 

Nodel 

0 1 

07 

02 

Node 2 

1 0 

00 

00 

Node 3 

1 0 

00 

00 


Table 3 17 Routing probabilities for class 1 customers of Example 3 7 


5 


































Node 1 

Node 2 

Node 3 

Node 1 

0 1 

07 

02 

Node 2 

1 0 

00 

00 

Node 3 

1 0 

00 

00 


Table 3 18 Routing probabilities for class2 customers of Example 3 7 
Output parameters for Example 3 7 from QNAT 



Node 1 

Node 2 

Node 3 

EN, 

0 886640 

0 886640 

4 22 

672 

classes 

Class 1 

Class 2 

Class 1 

Class 2 

Class 1 

Class 2 

EN.w 

044332 

0 44332 

0 44332 

0 44332 

2 1 13360 

2 11336 

Efy,k 

0 02310 

02310 

0 3300 

0 3300^ 

0 93800 

0 93800 

c ) 

8 67553 

8 67553 

6 07287 

6 07287 

1 735107 

1 735107 


1 0000 

1 0000 

0 7000 

0 7000 

0 2000 

0 2000 


Table 3 19 Output parameters of Example 3 7 


3 2 3 Mixed Net^\orks 

Mixed networks are those in which some classes are open and some are closed with 
respect to the system The open classes in the network are characterised b\ the total 
arrival rate of that class and the closed classes are characterised bv the population sector 
(population of each closed class) The algorithm is taken from the in in.rkii el „l In 
this algorithm the utilizations of the open classes at each node are computed first 
followed by the nett utilization of each node due to the open classes The pertormance 
parameters of the closed classes are computed by eliminating the open classes and 
converting the miked model into a closed model with inflated sen ice demands for the 

closed classes The inflation factor used IS 1 t/„o, (1 utilisation due to open classes at 

node i) and Ui (o) is given by 

U, (01 = L utilisation at node i of class c 


Graham and Kenneth C Se\cik Qiiaiut!att\e Swtcni 


' Edward D Lazowska JohnZahorjan G Scott 
Performance Prentice Hail Inc 1984 
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This IS the percentage of time the server is not used by open classes and is the time 
available for the closed classes to get service from the server in the queue The 
performance of the closed classes is analysed using the standard solution techniques 
discussed earlier with the inflated service times for the closed class customers The 
performance of the open classes are computed by considering the mean queue length of 
the closed classes at each node in the network The GUI permits the use of zero open 
class with non zero closed classes or vice versa 

Assumptions 

• All the queues in the network have infinite buffer capacity 

• The service discipline for all classes are FCFS at all the nodes in the network 

• There can be creation or combination of jobs for the open classes 

• All the classes are independent 

• There must be an external arnxal for each of the open classes These protes es 
must be Poisson 

• The ser\ ice time distribution for all the classes in the network are exponential 

• The routing chain for the jobs in the network are static 

Input parameters 

R total number of classes in the mixed model 
O total number of open classes in the mixed network 
C total number of closed classes in the mixed network 
R = 0 + C 

N number of nodes in the mixed model 

number of servers at node i in the mixed network for i = 1 2 
C population vector of the closed classes 
= { Cl C Cc } 

[QklvxA' qk(ij) routing probabilities for f;=l 2 J^&k=\2 R 


Open classes 

Xo, , mean arrival rate of class c at node i for c=l 2 O &. /=! 2 



T,c mean service time distribution of class c at node i cs. {0 ] 


forc=12 O & j=l 2 ^ 

Yic multiplication factor of class c at node i for i =1 2 &.c=\ 2 

Closed classes 

T,c mean service time distribution of class c at node i c e { C ) 
forc=l 2 C & t=l 2 ^ 


Output parameters 

H/Viopen mean open class jobs at node i 
ENt losed mean closed class jobs at node i 
Open classes 

X,, nett arrival of class c at the node i for /=1 2 jV <S. t-1 2 


jBNiopcn niean queue length of class c at the node i tor i- 1 2 


c-1 2 


mean delay of class c at the node i for /=! 2 N 


c=l 


0 


dep.open i departure rate of class c at the node / for 1 2 N 

c=l 2 O 

V.open Visit counts of class c at the node i for /= 1 2 N 


c=i 1 


O 


t/,, utilisation of class c at the node i for i=l 2 jV <S. c-1 2 
Closed classes 

losed c mean queue length of class c at node I for /=! 2 N 

c=l 2 C 

EW, losed mean delay of class c at node z for i=l 2 &c—\ 2 

XidosedL throughput of class c at node i for i=l 2 JV S. c-1 2 

Vi losed L 'Visit ratios of class c at node z for i=l 2 JV (S. c=l 2 



Algorithm 3 7 Mixed network solution technique 


R total number of classes 
O total number of open classes 
C total number of closed classes 
N number of nodes in the network 
{ 0 } set of open classes 
{ C } set of closed classes 

1) Solving the equation below we get the nett arrivals to the nodes of open classes 




^ik = k + X 0 1) 

/-.I 

for 7=1 2 & Jt=l 2 

0 

ikik 




f/.k= 


for 7= I 2 Jf &k=\ 2 

0 

771, 




0 




Ui (01 = X 

t=\ 

N 

for 7=1 2 A^&;:=12 

0 

dopiopcn k ~ ^ik Yik ( ^ 

1 

7=1 

qK(,j, ) for 7=1 2 N&,k=\2 

V 

throughput(A) = X 

X.k 

forA=12 0 



1=1 


2 ) 


V, for>=l2 .VSt=I2 0 

throughput(^) 

Inflate the service time of the closed classes b\ a factor of 1 Y,{o\ 


t.k = 


for 1=12 &. k=l 2 C & k E [C] 


1 Ut[0] 

Now solve the closed network with the algorithm 3 6 (multiple class closed 
network MVA ) 


' Edward D Lazowska John Zahorjan G Scott Graham and Kenneth C Sevcik QiamiiUitn e S^Mcm 
Performance Prentice Hall Inc 1984 
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( Algonthm 3 7 continued ) 


1 3) The waiting time and the queue length for the open classes are 
^lopcn k ( 1 +Qlength,jc) ) 

^.openk= fori=12 J^&h=\2 0&ke{0] 


1 [/, 


I (0) 


where 


Qlength, lo = EN, 

losed for 1=1 2 
c 

= X ^^lUosedk for 1=1 2 &ke [C] 


k-\ 


■E/Vjopen k = throughput(k) By.open k for i=l 2 7^^ &A=I2 0 & ke {0} 

0 

£'A'’.open=X ^^openk fOr 1=1 2 M &ke[ 0 ] 


*-l 


Example 3 8 




05 

▼ 

05 ^ 

(1) -121, 

I— r h 


A 

^ P ■ = ^ T 

s 

Figure 3 7 Multiple classes in Mixed Network (2 classes) 

(one open class and one closed class) for Example 3 8 


Description of the network for Example 3 8 

There are two nodes and two classes (one open class and one closed class) in the network 
Each open class is characterised by the total arrival rate (sum of arrival rate at node 1 and 
node 2) and each closed class is characterised bv the total number of customers in the 
network (population of the closed class is 6) The service time distributions at all nodes 
for all classes are exponential The routing probabilities and the service times at each 
node for each class are given below 





Input parameters to Example 3 8 


Open class 



Nodes 1 

Node 2 

^1 k 

0 500 

0 500 

T,k 

0 500 

0 500 

Y.k 

1 0 

1 0 


Table 3 20 Input parameters of open class in Example 3 8 


Closed class 



Nodes 1 

Node 2 


0 040 

0 280 


Table 3 2 1 Input parameters of closed class in Example 3 8 


Routing probabilities 


Ok Ml 

Nod 

esl 

Node 2 

Classes 

Class 1 

Class 2 

Class 1 

Class 2 

Node 1 

0 1 

02 

05 

08 

Node 2 

05 

0 8 

0 1 

02 


Table 3 22 Routing probabilities for Example 3 8 


Output parameters for Example 3 8 from QNAT . 
Open class 



Node 1 

Node 2 

f-N open 

1 944430 

11 38890 

k 

0 50 

0 50 

X,k 

1 125 

1 125 

f^iopen k 

1 944430 

11 38890 

£3V,oDen k 

1 944430 

11 38890 

dePiooen k 

0 50 

0 50 

Vcnenk 

1 12o 

I 125 

1 k 

0 6250 

0 6250 































Closed class 



Node 1 

Node 2 

■Enclosed 

0 1666582 

5 833340 

EA/itlosed k 

0 1666582 

5 833340 

EW iclosed k 

0 1777234 

3 60893 

^1 loscd k 

1 339276 

1 339276 

Vitlosed k 

1 0000 

1 0000 


Table 3 24 Output parameters of closed class in Example 3 8 




Chapter 4 


Summary and Future work 

4 1 Summary 

In this thesis we have developed a Microsoft Windows based software package for 
queueing network analysis called QNAT In this package some of the widely used 
solution techniques are used to solve a specified queueing network This package can be 
used to analyse a network of queues with queues having finite buffer capacity or infinite 
buffer capacity For infinite buffer capacit\ queues the jobs can be of a single class or oi 
multiple classes To solve a single class open network we have used GI/G/m algorithm 
and to solve a closed network we have used the state dependent MVA algorithm QNAT 
also allows us to include fork join queues that mav or may not have synchronising queue 
alter the service queue in the siblings We have developed reduction techniques to handle 
fork join queues without synchronising queues for open and closed networks QNAT 
does not support fork join queues with synchronising queues in an open network 

There are a few software packages that ha\e alreadv been de\ eloped to analvse queuein, 
networks We overcome the major drawback in these packages bv providing a use 
fnendly GUI based front end to define the network and to view the results of the solution 
Also the use of Mathematica™ to carrv out the computational tasks considerabK 
simplifies the developement process and enables us to use mathematically invohea 
algorithms Another salient feature of QNAT is that we have included many types or 
queueing networks in one package All the other packages that have been developea 
before restrict themselves to a small class of queueing networks 
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4 2 Future Work 

There are still many models that are not yet supported in QNAT Some of these are 
Norton reduction change of job classes in the case of multiple classes permit multiple 
classes in fork join queues permit multiple servers in the fork join model nodes with 
finite capacity m the fork join queues and nodes with finite capacity in the networks with 
multiple classes jobs Apart form the mathematical analysis of a queueing network there 
IS a need to provide simulation capability A simulation option when provided in QNAT 
can be used to define and execute a simulation model of the system under studys This 
simulation capability will also enable us to check the accuracy of the queueing network 
model We would also like to include discrete time queues and queueing networks into 
QNAT This IS expected to be particularly important in modeling communication 
networks based on ATM technology Allowing a combination of continuous time queues 
and discrete time queues in a network would be an interesting challenge 
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